Realization of a distributed Bragg reflector for propagating guided matter waves 
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We report on the experimental realization of a Bragg reflector for guided matter waves. A Bose-Einstein 
condensate with controlled velocity distribution impinges onto an attractive optical lattice of finite length and 
directly probes its band structure. We study the dynamics of the scattering by this potential and compare the 
results with simple one-dimensional models. We emphasize the importance of taking into account the gaussian 
envelope of the optical lattice which gives rise to Bragg cavity effects. Our results are a further step towards 
integrated atom optics setups for quasi-cw matter waves. 



The interaction of light with structures having a periodic 
refractive index profile is ubiquitous in photonics. Applica- 
tions range from simple antireflection coatings to the fabri- 
cation of dielectric mirrors with ultra-high reflectivities, used 
for instance in high-finesse cavities, and to semi-conductor 
laser technology with the example of VCSELs, and DFB or 
DBR lasers. In the field of guided optics, fiber Bragg gratings 
are essential components for the telecommunication industry, 
as well as for the realization of outcoupling mirrors in high- 
power fiber lasers. Photonic crystal devices also have a huge 
range of applications [1]. 

In matterwave optics and interferometry, interactions of 
free-space propagating beams or trapped clouds with periodic 
structures or potentials have been extensively investigated and 
are commonly used as mirrors and beamsplitters [2] . In this 
paper, we demonstrate, following the proposals of Refs. [3- 
7], a Bragg reflector for manipulating a guided Bose-Einstein 
condensate (BEG) propagating in an optical waveguide, i.e. 
the exact atom-optics counterpart of a photonic fiber Bragg 
grating. We study the dynamics and the transmission of a 
probe wavepacket as a function of its incident velocity and 
of the depth of the optical lattice. As we will develop later on, 
this quasi ID configuration exemplifies clearly two textbook 
features of quantum mechanics: quantum reflection [2, 9, 10] 
and band theory [11, 12]. This article is organized as fol- 
lows. We first present a simple model to gain some physical 
insight into the Bragg reflection of a matter wavepacket by 
a finite-length lattice having a gaussian envelope. Then we 
describe our experimental implementation and show quantita- 
tive agreement between the data and our model. Finally, we 
discuss numerical simulations that give access to unresolved 
details in the experiment. 

We consider a BEG with given mean velocity v and disper- 
sion Av propagating in a horizontal waveguide defining the 
X-axis. At some distance, two intersecting laser beams inter- 
fere and create an attractive quasi-periodic potential of finite 
length, with lattice spacing d (see Fig. 1). 

The potential experienced by the atoms is modeled by: 
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whose depth Uq is proportional to the power of the laser 
beams. The period d naturally introduces typical scales in 



velocity vr = h/{md) and energy E^i = mv^/2. 

We are interested in a wavepacket with finite size and ve- 
locity dispersion impinging on a finite-length lattice. Let us 
consider first the textbook case of an incident plane wave and 
a square envelope (see e.g. [13] for an analytical treatment of 
the problem). Figure 2(a) shows the transmission coefficient 
for a lattice of N = 800 sites, calculated by solving numeri- 
cally the corresponding stationary Schrodinger equation. For 
a given velocity v, one observes that the transmission coef- 
ficient essentially switches between and 1 as a function of 
the depth of the lattice Uq. It can be interpreted as follows. 
At the entrance, the incoming state of energy Ei = mv'^ /2 
is projected onto the eigenstates of the lattice. The associ- 
ated eigenenergies distribute into energy bands [11] whose 
position can be expressed in terms of the characteristic func- 
tions of the Mathieu equation [1, 15] as depicted by the white 
dashed lines in Fig. 2(a). Reflection occurs if Ei lies in the gap 
between two allowed energy bands. Due to the finite length 
of the lattice the energy bands are not strictly continuous and 
resolve into N separate states for vanishing incident veloc- 
ity [15]. Undersampling of the image gives rise to the 'foamy' 
aspect of the low- velocity side of transmission bands. Obvi- 
ously, the reflection by an attractive potential is a purely quan- 
tum effect, with no classical counterpart. 

A second interesting feature appears in the limit of a van- 
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FIG. 1 : (color online) Schematic view of the experimental setup (not 
to scale). 
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FIG. 2: (color online) Transmission coefficient of the lattice as a 
function of the incident velocity v and of the lattice depth Uq. (a): 
Square-envelope lattice with 800 sites. The white dashed lines are 
obtained from the Mathieu characteristic functions; white arrows 
show the velocities for which the Bragg condition is fulfilled (see 
text), (b): Lattice with a gaussian envelope (1/e^ radius w ^ 230d). 
The vertical shaded stripe corresponds to the relevant velocity com- 
ponents in the wavepacket used for the measurements shown in 
Figs. 3 and 4. The insets on top of panels (a) and (b) illustrate the 
shape of the lattice enveloppe, but the number of sites is reduced to 
N = 20 for clarity. 



ishing potential depth Uq. One still observes reflection but it 
occurs only for incident velocities of the form v = nvji/2 
where n is an integer (see the white arrows on Fig. 2(a)). This 
corresponds to Bragg's condition 2dsmQ = nX [12], where 
9 = 7r/2 for retroreflection, and A = h/{mv) is the incident 
de Broglie wavelength of the atoms: the reflection amplitude 
at each lattice site is small but constructive interference be- 
tween all the reflected waves results in a macroscopic reflected 
wave building up. For the range of parameters of Fig. 2(a), 
quantum reflection by a single lattice well occurs only for ve- 
locities that are very small as compared to vr [15] and thus 
cannot explain the observed features. 

Figure 2(b) shows the result of the same calculation, but 
now for the experimentally relevant case of a lattice having 
a smooth gaussian envelope. One clearly observes a drastic 
change in the dependence of the transmission coefficient: for 
a fixed velocity v, the transmission is essentially equal to one 
only below a critical value Uq ^ of the lattice depth, and then 
goes to zero. That critical value corresponds to the smallest 
one at which total reflection would occur for the square enve- 
lope lattice. No resurgence of the transmission is observed if 
Uq is further increased which gives a 'sawtooth' appearance 
to the boundary between reflection and transmission. 

This can be understood as follows. We are in the slowly 
varying envelope limit sls w :^ d. The amplitude of the lat- 
tice does not change appreciably over a few sites, and thus 
one can consider that, locally, the matterwave interacts with 
a constant-amplitude lattice. When Uq = Uq the reflection 
condition is met at the center of the lattice i.e. at x = 0. Then, 
when Uo increases, there are some locations ixrefi, on both 
sides of the center, for which /7(±Xrefl) = Uq ^ and where re- 
flection occurs. This constitutes a Fabry-Perot resonator made 
of two Bragg mirrors, analog to VCSELs for example. As in 



optics, transmission exhibits sharp resonances which gives the 
same foamy aspect as in Fig. 2(a). 

We now come to the experimental realization. Our tech- 
nique to produce all-optical BECs has been described in de- 
tails elsewhere [16]; in what follows we thus simply recall 
the major steps. We produce an almost pure ^^Rb condensate 
containing typically 5 x 10^ atoms by forced evaporation over 
4 s in a crossed optical dipole trap. It is made of two inter- 
secting beams with a wavelength of 1070 nm. A horizontal 
one, with a waist of 50 /im, to be used later as a guide for 
the BEC, defines the x direction. The second, the 'dimple' 
beam, of waist 150 /im, propagates along the f + i direction, 
z being the vertical (Fig. 1). Spin distillation using a magnetic 
field gradient during evaporation [16] results in the BEC be- 
ing prepared in the state |F = 1, tuf = 0). We then decrease 
adiabatically the power in the dimple beam by a factor ~ 20 
over 80 ms, thus barely keeping a longitudinal confinement 
for the BEC, before switching it off abruptly to outcouple a 
wavepacket in the horizontal guide. In this way, we produce a 
wavepacket with a small longitudinal velocity dispersion. To 
set the wavepacket in motion, we then switch on a coil, coax- 
ial with the guide, that produces an inhomogeneous magnetic 
field. Through the quadratic Zeeman effect the wavepacket is 
accelerated for typically 15 ms to a final mean velocity v be- 
tween 2 and 15 mm/s. The residual acceleration of the packet 
due to stray fields and beams curvature is negligible (we mea- 
sure an upper bound of 10 mm/s^). 

Centered 250 /im downstream from the dimple location xq, 
the optical lattice is produced at the intersection of two beams 
with a wavelength Al = 840 nm (red-detuned with respect to 
the ^^Rb Dl and D2 lines) and a waist w = 110 /im, linearly 
polarized along i, crossing at an angle c^ 81°. The lattice 
detuning is large enough so that spontaneous emission does 
not play any role on our experimental timescales. The result- 
ing lattice spacing is d = Al/[2 sm{6/2)] ~ 650 nm. In a set 
of preliminary experiments we calibrate the potential depth Uq 
using Kapitza-Dirac diffraction [17, 18]. A BEC is created at 
the position of the lattice and exposed to the lattice potential 
for a short time tkd, typically a few tens of microseconds. 
The diffraction pattern of the BEC after time of flight as a 
function of tkd is then compared to numerical simulations of 
the process. A typical 35 mW per beam results in Uq up to 

15^R. 

After being launched as described above, the wave packet 
propagates in the horizontal guide for an adjustable time tprop- 
Then all the lasers are switched off abruptly and the cloud is 
imaged by absorption after a 10 ms time of flight. This gives 
access to the spatial density distribution n(x, t) = |V^(x, t)p 
of the wavepacket with a resolution of about 10 /im limited by 
the numerical aperture of our collection lens. 

In a first set of experiment, the propagation time tp^op = 
100 ms is sufficiently long so that the interaction with the lat- 
tice is completed. We measured in a separate experiment a 
mean velocity is ^ c^ 11 mm/s c^ 1.6 vr and a r.m.s veloc- 
ity spread Av c^ 1.3 mm/s c^ 0.2 vr corresponding to the 
shaded region of Fig. 2(b). For each lattice depth Uq, an aver- 
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FIG. 3: (color online) (a) Measured density distribution of the 
wavepacket (of initial mean velocity v = 11 mm/s) after a prop- 
agation time tprop = 100 ms, for different lattice depths Uq. Each 
horizontal line is the average of typically eight absorption images in- 
tegrated along the y direction. The vertical white dashed line shows 
the position of the center of the lattice, (b) Results of the simula- 
tion without any adjustable parameters. The finite resolution of the 
imaging system (^10 /xm) is included. The dotted line is the ex- 
pected position of the center of the wavepacket according to classical 
dynamics (see text). 
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age image is generated from eight individual runs and then in- 
tegrated along the transverse direction y. Figure 3(a) is a stack 
of 55 such profiles. For sake of comparison, Figure 3(b) is the 
result of a numerical simulation of the wavepacket dynam- 
ics using the one-dimensional Schrodinger equation solved 
by the split-Fourier method; the initial condition is a gaus- 
sian wavepacket with the experimentally measured momen- 
tum and position dispersions [19]. There is no adjustable 
parameter and the overall agreement with experimental data 
means that our simple ID model captures most of the physics 
involved. 

Let us concentrate first on the transmitted part of the 
wavepacket (x > 0). If there were no lattice, the propaga- 
tion time tprop is long enough so that the initial size of the 
wavepacket is negligible with respect to its size after propa- 
gation. The spatial distribution of the wavepacket would then 
be a direct mapping of its initial velocity distribution f{v)\ 

n{x, tprop) OC f[{x- Xo)/tprop]. 

Figure 3(b) compares then directly with Fig. 2(b) whose 
horizontal axis is properly scaled. In the background of the 
shaded area of Fig. 2(b) representing the wavepacket one can 
see the transmitted and reflected components. In the presence 
of the lattice, the reflected part propagates backwards and is 
located, for the propagation time chosen here, at a symmetri- 
cal position. This explains why the transmitted and reflected 
wavepackets appear like complementary mirrored image of 
each other. The sawtooth-like boundary, reminiscent of the 
transmission diagram, is a fingerprint of the band structure in- 
side the lattice. However, the effect of the lattice potential is 
not limited to the one of the sinusoidal component, respon- 
sible for the Bragg reflection described above. The spatially 
averaged attractive potential also accelerates the wavepacket. 
The white dotted line on Fig. 3(b) shows the final position 
of a classical particle starting with velocity v from position 
xq and propagating for a time tprop, taking into account its 
acceleration by the spatially averaged lattice potential. The 



FIG. 4: (color online) Time sequence showing the scattering of a 
wavepacket with mean velocity i; ^ 1 1 mm/s on the optical lattice 
for Uq ^ ll^R. The white dashed lines in panels (a) and (b) show 
the position of the center of the lattice. The time interval between 
successive images is 1 ms. (a): Experimental data, (b): Simulation, 
taking into account the finite imaging resolution as well as the time 
of flight (TOF) period, (c): Same as (b) but without TOF nor reduced 
resolution; the color scale is nonlinear in order to enhance contrast. 



fair agreement with the data indicates that the slight curvature 
in the position of the wavepacket as a function of Uq simply 
arises from this classical effect. 

Beyond studying the asymptotic scattering states, it is also 
possible to visualise the dynamics of the interaction by vary- 
ing tprop. Figure 4(a) displays such a time sequence that 
fairly compares to the numerical simulation depicted in the 
same conditions in Figure 4(b). One clearly observes the 
spreading of the incident wavepacket over the whole lattice 
for 30 ms < tprop ^ 45 ms and its subsequent splitting into 
a reflected and a transmitted one. Unfortunately, the details 
of the inner dynamics are washed out by the free expansion 
of the wavepacket during the time-of-flight sequence and the 
finite resolution of the imaging system. 

Numerical simulations, properly checked against the pre- 
vious experimental results, are useful here. In Figure 4(c) 
we have deliberately suppressed the time-of-flight period and 
enhanced the optical resolution and the contrast with respect 
to Fig. 4(b): one then clearly observes multiple reflections 
of some components of the wavepacket at symmetric posi- 
tions ixrefl, with decreasing amplitude at each bounce. This 
'cavity-ring-down' behavior explains the formation of struc- 
tures in the transmitted and reflected wavepackets as observed 
in Fig. 4(a) and especially visible as a parallel lower stripe in 
the transmitted wavepacket for 50 ms < tprop ^ 75 ms. 
However, experimentally, observing several bounces is not 
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the guide, in order to enhance the effects of nonUnearities. 
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FIG. 5: (a) Calculated density and (b) velocity distribution of a 
wavepacket of initial mean velocity ^; = 11 mm/s after a propa- 
gation time tprop = 100 ms with full resolution. The lattice depth is 



possible here due to the small number of atoms involved. 

In the same way, on images such as Fig. 3 the reflected 
wavepacket appears to be relatively smooth. It is actually not 
the case as can be seen on simulations with full resolution 
(Fig. 5). The lattice acts as a matterwave interference filter 
with very narrow features due to the high number of lattice 
sites (foamy zones of Fig. 2). 

Until now we have used a simple one-dimensional descrip- 
tion of the system. However the system is actually far from 
being one-dimensional, since the transverse quantum of en- 
ergy huj± 0:^ /i X 90 Hz is much smaller than the typical 
longitudinal energy scales, typically by two orders of mag- 
nitude. Our simple one-dimensional model agrees well with 
the experimental results as shown above because couplings 
between longitudinal and transverse degrees of freedom are 
weak (they are due only to experimental imperfections such as 
misalignments of the lattice beams with respect to the guide 
for instance); some transverse excitations can nevertheless be 
observed on our data (see e.g. the long wavelength dipole os- 
cillations on Fig. 4(a), especially for tprop ^ 30 ms). Stronger 
couplings would be expected to alter significantly the scatter- 
ing properties of the structure [20, 21]. 

In conclusion, we have studied in details the scattering of 
a guided matterwave by a finite length optical lattice in the 
slowly varying envelope limit. The experiments can be in- 
terpreted in the framework of a local band structure and the 
whole lattice can be seen as a Bragg reflector/cavity. 

Major improvements are expected with the use of high nu- 
merical aperture optics [22, 23]. Reducing drastically the 
length w of the lattice and thus generating a structure consist- 
ing of only a few sites, possibly with a shaped envelope, one 
could tailor almost arbitrarily the matterwave filter response. 
Moreover, if a controlled frequency offset between the two 
interfering beams is introduced, one generates a moving lat- 
tice. The transmission band of the filter could then be adjusted 
at will. Such setups would prove useful in measuring for in- 
stance the coherence length [24] of guided atom lasers [25- 
27]. In a different direction, it would be appealing to study 
the effect of interatomic interactions [28] on the propagation 
of the wavepacket, with the possible appearance of soliton 
trains [29], or atom-blockade effects [30]. This regime could 
be reached by using much higher transverse frequencies for 
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Supplementary material 



In this supplementary material, we provide additional details concerning several points discussed in the main text about the 
physics of Bragg reflection by a finite-length, square-envelope lattice. 



TRANSMISSION BY A SEMI-INFINITE, 

SQUARE-ENVELOPE LATTICE AND MATHIEU 

CHARACTERISTIC FUNCTIONS 

Here we detail how one can obtain analytically, using the 
theory of the Mathieu equation, the regions in parameter space 
{v^Uo) corresponding to total reflection by a semi-infinite lat- 
tice. 

The time-independent Schrodinger equation for a parti- 
cle of mass m and energy mv'^/2 evolving in the potential 
— Uo sin^(7rx/(i) reads: 



— Jl^ d'^l/j . 2 fTTX 

2m dx^ V d 
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^{x). (2) 



Introducing the dimensionless parameter x = irx/d, Eqn. (2) 
can be rewritten as 



dx^ 



[a — 2gcos(2x)] '0 = 0, 



(3) 



which is precisely the definition of the Mathieu equation [1], 
with a and q defined by 



a = 4{v/vnf^2Uo/E^, 
q = Uo/E^. 



(4) 



The theory of the Mathieu equation (3) shows that, de- 
pending on the parameters (g, a), either the solutions remain 
bounded whatever the initial conditions (stable regions of the 
(g, a) plane), or there is at least one solution that grows ex- 
ponentially when X -^ oo [1]. The stable regions in the 
(g, a) plane are shown in red in Fig. SI (a). The equation of 
the boundaries between stable and unstable regions are given 
by the so-called Mathieu characteristic functions an{q) and 
^n(^). where n G N. 
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Fig. SI: (a): Stability diagram of the Mathieu equation in the (^, a) 
plane. The stable regions correspond to the red areas, (b) The same 
diagram as a function of the parameters {v,Uo), see Eq. (4). 



For the problem of a quantum particle impinging with ve- 
locity V on a semi-infinite lattice of depth Uq, one easily finds, 
using the mapping (4), the regions in the (v^Uq) parameter 
space that correspond to unstable solutions of the Mathieu 
equation. They correspond to situations in which the energy of 
the incident particle lies in a forbidden band, implying a total 
reflection by the lattice. These are shown in blue on Fig. Sl(b) 
and agree well with the regions of full reflection by a finite- 
length lattice shown in the main text on Fig. 2(a). 

However, for parameters corresponding to full transmission 
(in red on Fig. Sl(b)), and especially at low velocities, one 
clearly sees on Fig. 2(a) of the main text that for a finite-length 
lattice, significant reflections are observed in the regions with 
a 'foamy' aspect. These arise from the finite number of lattice 
sites, as we shall see in the section below. 



INFLUENCE OF THE NUMBER OF LATTICE SITES FOR A 
SQUARE-ENVELOPE LATTICE 

The above analytical results are valid for a semi-infinite 
lattice. For the more experimentally relevant case of a 
finite-length lattice, we solve numerically the stationary 
Schrodinger equation. Figure S2 shows, for increasing num- 
bers of sites N of the lattice, the calculated transmission as a 
function of (v^Uq). 

Figure S2(a) corresponds to the transmission of a single si- 
nusoidal site. One clearly observes quantum reflection at low 
velocities (typically, below 0.3 vr for the range of depths Uq 
studied here), except for discrete values of U^ that correspond 
to the appearance of a new bound state in the attractive well: 
in that case instead, full transmission occurs. Note that in or- 
der to observe such a quantum reflection by a single lattice site 
with the velocities of 1.6 vr used in the experiments described 
in the main text, one would need to have a much deeper well 
(typically several 10^£^r). 

When increasing the number N of sites, two main effects 
arise (Fig. S2(b)-(e)): 

• At low velocities, the discrete values of Uq correspond- 
ing to transmission split into exactly N components, 
due to the coupling between adjacent sites [2]. When 
the number of sites becomes large, these transmission 
resonance spread all over the allowed energy bands 
given by the Mathieu equation as can be seen by com- 
paring Figs. S2(e) and Sl(b). The width of these reso- 
nances thus scales as 1/N with the number TV of sites. 
For our experimental parameters (TV ~ 800), they are 
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Fig. S2: Transmission coefficient of a square-envelope lattice of (a)-(e): N = 1, 2, 3, 5, 25 sites, (f): Magnification of the region in the white 
rectangle in (e), showing clearly that the allowed "band" is actually made of N narrow resonances. 



extremely narrow and can be observed only in the nu- 
merical simulations (see Fig. 5 of the main text). 

• The regions in the {v^Uq) plane in which total reflec- 
tion occurs extend towards higher velocities, approach- 
ing the ones obtained in the limit A^ ^ oo in the section 
above: in that case, the reflection originates from Bragg 
constructive interference of the waves scattered by mul- 
tiple sites, and not from quantum reflection on the steep 
slope of a single attractive potential. 
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